Dangling-bond spin relaxation and magnetic 1// noise from the 
amorphous-semiconductor/oxide interface: Theory 
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We propose a model for magnetic noise based on spin-flips (not electron-trapping) of paramagnetic 
dangling-bonds at the amorphous-semiconductor/oxide interface. A wide distribution of spin-flip 
times is derived from the single-phonon cross-relaxation mechanism for a dangling-bond interacting 
with the tunneling two-level systems of the amorphous interface. The temperature and frequency 
dependence is sensitive to three energy scales: The dangling-bond spin Zeeman energy (S), as well 
as the minimum (-Emin) and maximum (-Emax) values for the energy splittings of the tunneling two- 
level systems. At the highest temperatures, ksT 3> Max (5, i?max), the noise spectral density is 
independent of temperature and has a 1// frequency dependence. At intermediate temperatures, 
ksT ^ 5 and -Emin ^ ksT ^ Emax, the noise is proportional to a power law in temperature 
and possesses a 1/f^ spectral density, with p = 1.2 — 1.5. At the lowest temperatures, ksT <^ 5, 
or fc^r <C Emin, the magnetic noise is exponentially suppressed. We compare and fit our model 
parameters to a recent experiment probing spin coherence of antimony donors implanted in nuclear- 
spin-free silicon [T. Schenkel et ai, Appl. Phys. Lett. 88, 112101 (2006)], and conclude that a 
dangling-bond area density of the order of 10^^ cm~'^ is consistent with the data. This enables the 
prediction of single spin qubit coherence times as a function of the distance from the interface and 
the dangling-bond area density in a real device structure. We apply our theory to calculations of 
magnetic flux noise affecting SQUID devices due to their Si/Si02 substrate. Our explicit estimates 
of flux noise in SQUIDs lead to a noise spectral density of the order of 10~^^$o(Hz)^^ at / = 1 Hz. 
This value might explain the origin of flux noise in some SQUID devices. Finally, we consider the 
suppression of these effects using surface passivation with hydrogen, and the residual nuclear-spin 
noise resulting from a perfect silicon-hydride surface. 

PACS numbers: 05.40.Ca; 61.43.-j; 76.30.-v; 85.25.Dq. 



I. INTRODUCTION 

Our physical understanding of spin relaxation in semi- 
conductors plays a crucial role in the current develop- 
ment of spin-based electronics ""^ and spin-based quantum 
computation.^ One question that received little or no at- 
tention so far is related to magnetic noise in semicon- 
ductor devices and nanostructures. Magnetic noise from 
impurities and other defects at the interface may be the 
dominant source of spin phase relaxation (decoherence) 
for implanted donor electrons'' or nuclear spins^ in iso- 
topically purified silicon. Moreover, because Si/Si02 and 
other amorphous oxide interfaces are used as the sub- 
strate for sensitive SQUID magnetometers,^"^ the spin 
relaxation of magnetic impurities at the substrate might 
explain the observed magnetic flux noise in these devices. 

One universal characteristic of silicon devices is the 
presence of an insulating interface, usually an oxide, sep- 
arating the metallic gate from the semiconductor. It is 
known for a long time that these interfaces are rich in 
dangling-bond type defects (also denoted "Ph centers") 
which can be detected using spin resonance techniques. 
These studies have established a wide distribution of 
dangling-bond (DB) energy levels, spanning almost the 
whole semiconductor energy gap, with each DB charac- 
terized by a large on-site Coulomb energy U ^ 0.5 eV.®'^ 
When the dangling-bond (DB) energy level falls within 
ksT of the interface Fermi level, it acts as a trapping- 
center and leads to the well known 1// charge and cur- 



rent noise for interface conduction electrons. Neverthe- 
less at low temperatures the area density for trapping- 
center DBs is only a tiny fraction of the area density for 
paramagnetic DBs. For example, at T = 5 K this frac- 
tion is only ksT/U ^ 10^^ (Fig- !)• As a consequence, 
the magnetic noise due to paramagnetic DBs is at least 
a factor of U/ksT ^ 1 larger than magnetic noise gen- 
crated by electron trapping, provided the paramagnetic 
DBs have a non-zero spin-flip rate (Magnetic noise due 
to electron trapping is discussed in appendix A). 

The spin relaxation rate for dangling-bond type defects 
depends crucially on the non-crystalline nature of amor- 
phous compounds. ^^"^"^ However, a detailed theoretical 
study of the magnetic field and temperature dependence 
of this effect has not been done. In this article we present 
a general theory of dangling-bond spin-lattice relaxation 
in amorphous materials, and show that the noise created 
by the magnetic dipolar field of an ensemble of dangling- 
bonds has the 1// frequency dependence at high temper- 
atures. We fit our theory to a recent experiment probing 
spin coherence of antimony donors implanted in nuclear- 
spin-free silicon'^ in order to estimate our model param- 
eters. 

We exploit the important relationship between phase 
coherence of a localized "probe" spin (e.g. the implanted 
Sb spins in Ref. 3) and its environmental magnetic noise 
(Fig. 2). The coherence decay envelope of a "probe" spin 
measured by a class of pulse spin resonance sequences 
is directly related to a frequency integral over magnetic 
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FIG. 1: (Color online) Band diagram for a Si/Si02 interface. 
Dangling bonds with energy much larger than ep are empty; 
DBs with energy in the interval (ef — ksTjep + ksT) are 
trapping-centers for interface conduction electrons, responsi- 
ble for charge, current, and magnetic noise. DBs with energy 
in the interval (ep — U,eF — ksT) are singly occupied (param- 
agnetic), and hence contribute exclusively to magnetic noise. 
DBs with energy less than ep — U are doubly occupied and 
do not contribute to any kind of noise. 



noise times a filter function. This allows us to interpret 
pulse spin resonance experiments of localized spins as 
sensitive detectors of magnetic noise in nanostructures. 
The spin qubit phase coherence is a local probe of low 
frequency magnetic noise. The same ideas apply equally 
well to experiments probing the coherent dynamics of 
superconducting devices. ^^'^^ 

An important step towards this characterization was 
given recently, by the report of the the first measurements 
of spin echo decay in silicon implanted with an ultra-low 
dose of antimony donors (~ 10^^ cm~^).'^ Two samples 
were reported, 120 KeV and 400 KeV, with low and high 
implant energy respectively. The former leads to a donor 
distribution closer to the interface, see Table I. 

Table I provides experimental evidence that the sur- 
face leads to additional mechanisms for donor spin phase 
fluctuation and magnetic noise. These mechanisms seem 
to contribute exclusively to the phase coherence time 
{T2) but not to the spin-flip time (Ti) of the Sb donors, 
therefore the associated noise spectrum should be low 
frequency in nature (with a high frequency cut-off much 
smaller than the spin resonance frequency). 

Here we consider the mechanisms of magnetic noise 
that might be playing a role in these experiments. For 
a Si/Si02 interface we show that dangling-bond spin- 
flips play a dominant role. A dangling-bond (DB) is a 
paramagnetic defect usually associated with an oxygen 
vacancy in the Si/Si02 interface. These point defects are 



Sample 


Interface Peak depth [nm] 


Ti [ms] 


T2 [ms] 


120 KeV 


Si/Si02 


50 


15 ±2 


0.30 ±0.03 


120 KeV 


Si-H 


50 


16 ±2 


0.75 ±0.04 


400 KeV 


Si/Si02 


150 


16 ± 1 


1.5 ±0.1 


400 KeV 


Si-H 


150 


14 ± 1 


2.1 ±0.1 



TABLE I: Spin relaxation data'^ taken at 5.2 K for antimony 
donor electron spins implanted in isotopically purified silicon. 
Ti was measured using inversion recovery ESR, while T2 is the 
1/e decay of Hahn echo. For each sample, data was taken for 
the untreated oxidized surface (Si02) and for the passivated 
surface, treated with hydrofluoric acid in order to obtain a 
hydrogen terminated surface. The data clearly indicates that 
(1) donors close to the surface have lower spin coherence times 
T2 but the same spin-flip time Ti; (2) Surface passivation leads 
to a sizable increment in T2, but no change in T\. 



generically denoted "P;, centers" with chemical structure 
represented by Sis = Si-.^^^ There is yet no experimental 
or theoretical studies of spin relaxation times (Tf ^) for 
DBs at the Si/Si02 interface. Nevertheless a systematic 
study of DB spin relaxation in bulk amorphous silicon 
was carried out in the lOSO's.^^'^'^ The measured DB spin 
relaxation rate was found to increase as a power law on 
temperature, 1/Tf ^ cx T" with an anomalous exponent 
n = 2 — 4 dependent on the sample preparation method. 
At T = 5 K and B = 0.3 T the typical T^^ was in the 
range 0.1 — 1 ms.^^ 

At first it seems puzzling that the dangling-bond spin 
would relax in such a short time scale at the lowest tem- 
peratures. The typical Ti of localized electron spins 
in crystalline silicon (e.g. phosphorus donor impuri- 
ties) is almost a thousand seconds in the same regime. 
This happens due to the weak spin-orbit coupling in 
bulk crystalline silicon. However, dangling-bonds in non- 
crystalline silicon are coupled to unstable structural de- 
fects, and this fact seems to explain their short Ti.^^'^^ 
These structural defects behave as tunneling two level 
systems strongly coupled to lattice vibrations (phonons). 
Each time a tunneling two level system (TTLS) under- 
goes a phonon-induced transition, the DB spin feels a 
sudden shift in its local spin-orbit interaction, which 
may be quite large because the TTLS is associated with 
a local reordering of the atomic positions of the non- 
crystalline material. As a consequence, the DB spin may 
flip each time the TTLS switches. Remarkably, this 
cross-relaxation process remains effective even at zero 
magnetic field because it does not involve a Kramers con- 
jugate pair (in contrast to spin-flips without a simulta- 
neous TTLS switch). 

We develop this theory further in order to incorporate 
the exponentially wide TTLS parameter distribution typ- 
ical of amorphous materials. As a result, we flnd that the 
magnetization of an initially polarized ensemble of DB 
spins will undergo non-exponential relaxation in time. 
Our theory of dangling-bond spin-lattice relaxation and 
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magnetic noise is based on an effective Hamiltonian ap- 
proach, allowing us to draw generic conclusions about the 
frequency, temperature, and magnetic field dependence 
of spin- noise in a variety of amorphous materials. For 
example, our results apply equally well to the magnetic 
noise produced by E' centers in bulk Si02, another well 
studied dangling-bond. Other materials of relevance to 
our work are the bulk AI2O3 (sapphire), and AI/AI2O3 
and Si/Si3N4 interfaces, whose paramagnetic dangling- 
bonds/magnetic impurities are yet to be characterized 
experimentally. 

Our results are of particular importance to magnetic 
flux noise in SQUID devices, whose microscopic origin is 
a longstanding puzzle (for a review see section IV-G of 
Ref. 18). In section VII we apply our results to calcula- 
tions of flux noise due to DBs within the area enclosed by 
the SQUID loop, and show that this contribution might 
explain some of the available flux noise measurements. 

It is possible to considerably reduce the dangling- 
bond area density using a surface passivation technique. 
For example, the application of hydrofluoric acid to the 
Si/Si02 surface removes dangling-bonds by covering the 
surface with a monolayer of hydrogen atoms. Recently, 
Kane and collaborators fabricated a field-effect- transistor 
using a passivated Si(lll)H surface, and demonstrated 
record high electron mobility. Nevertheless, the large 
density of hydrogen nuclear spins might be an impor- 
tant source of magnetic noise. The nuclear spins are 
constantly fluctuating due to their mutual dipolar cou- 
pling. In section VIII we consider calculations of mag- 
netic noise due to a hydrogen terminated Si(100)H sur- 
face. We use the same theory previously developed for 
Hahn echo decay of a phosphorus impurity in bulk doped 
natural silicon. We show that the Hahn echo decay 
in a Si(100)H surface has many peculiarities, including a 
special crystal orientation dependence for the donor T2 
times that may be used as the fingerprint for detecting 
this source of noise experimentally. 



II. RELATIONSHIP BETWEEN MAGNETIC 
NOISE AND PHASE RELAXATION IN PULSE 
SPIN RESONANCE EXPERIMENTS: 
ELECTRON SPIN AS A LOCAL PROBE OF 
MAGNETIC NOISE 

Consider the following model Hamiltonian for the in- 
teraction of a localized spin with a noisy environment, 



n = -^eBcj^ + fi{t) ■ cr. 



(1) 



Here a = {ax, cry,az) is the vector of Pauli matrices de- 
noting the state of the electron spin being probed by a 

pulse spin resonance experiment (henceforth called the 
donor spin - e.g. the Sb spins in Ref. 3), jgB is the spin 
Zeeman frequency in an applied external magnetic field 
B, and 7e = ge/(2TOec) is a gyromagnetic ratio for the 
electron spin [for a group V donor impurity such as P or 



Sb, 7e « 1.76 X 10^ (sG)~^ is close to the free electron 
value]. Note that Eq. (1) was divided by h so that en- 
ergy is measured in units of frequency. Each component 
of the vector i) = {fix^fly^flz) is an operator modeling the 
magnetic environment (the DB or other impurity spins) 
surrounding the donor spin. The simplest way to describe 
the time evolution of the spin's magnetization (cr) is the 
Bloch-Wangsness-Redfield approach, which assumes (cr) 
satisfies a first order differential equation in time. The 
decay rate for (cr^) is then given by 



1 



\ E [^«(+7eB) + 5,(-7eS) 



(2) 



q=x,y 



with the environmental noise spectrum defined by 

1 /"^ 

Sii^) = ^J_ e^^'iUm.mdt. (3) 

Note that the energy relaxation time Ti for the donor 
spin is determined by the noise at a; = ±7gi?, that is just 
a statement of energy conservation. Within the Bloch- 
Wangsness-Redfield theory the spin's transverse magne- 
tization ((cr_|_) = (cTj, +iay)/2) decays exponentially with 
the rate 



^ = ^+^4(0), 



(4) 



where wc added a * to emphasize this rate refers to 
a free induction decay (FID) experiment. The Bloch- 
Wangsness-Redfield approach leads to a simple exponen- 
tial time dependence for all spin observables. Actually 
this is not true in many cases of interest, including the 
case of a group V donor in bulk silicon where this approx- 
imation fails completely (for Si:P the observed Hahn echo 
decay fits well to e^"^ in many regimes). The prob- 
lem lies in the fact that the Bloch-Wangsness-Redfield 
theory is based on an infinite time limit approximation, 
that averages out finite frequency fluctuations. Note that 
T2* differs from Ti only via static noise, Sz{0) in Eq. (4). 
A large number of spin resonance sequences, most no- 
tably the Hahn echo are able to remove static noise com- 
pletely. 

We may develop a theory for spin decoherence that 
takes into account low frequency fluctuations in the 
semiclassical regime <C ksT, when Sz{—(^) = 
Q-f^/^BT g^^^^ « Sz{uj). The spin coherence envelope 
may be calculated in the pure dephasing limit (17a, = % = 
0), with the assumption that t)^ rjz is distributed ac- 
cording to Gaussian statistics. For derivations and dis- 
cussions on the applicability of this theory, we refer to 
Ref. 14. A similar method in the context of supercon- 
ducting qubits was proposed in Ref. 16. The final re- 
sult is a direct relationship between phase coherence and 
magnetic noise according to 



|(c7+(i))|=exp 



duj Sz{co)J^{t,w) 



(5) 
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FIG. 2: (Color online) How to detect low frequency magnetic 
noise using electron spin resonance. A low density of anti- 
mony (Sb) donor impurities is implanted in a Si/Si02 sam- 
ple using an ion gun, and the distribution of Sb donors is 
determined using secondary ion mass spectroscopy. Next, a 
Hahn echo decay experiment is performed on the Sb spins. 
The Hahn echo decay envelope is directly related to magnetic 
noise produced by e.g. dangling bonds at the interface, see 
Eq. (5). 

with !F{t, uj) a filter function that depends on the partic- 
ular pulse spin resonance sequence. For a free induction 
decay experiment {tt /2 — t — measure) we have 

.fFID(i,07) = -— ^— ^, (6) 
^ (cj/2) 

while for the Hahn echo (7r/2 — r — tt — r — measure) the 
filter function becomes 

1 sin^ (a>r/2) 

•?^Halm(2r,tj) = - o . (7) 

^ (ijj/4) 

Note that in the limit t — > cx) Eq. (6) becomes ■n5{Lo)t, 
recovering the Bloch-Wangsness-Redfield result Eq. (4). 
The Hahn echo filter function satisfies .FHahn(2r, 0) = 0, 
showing that it filters out terms proportional to 5*2(0) 
in spin evolution. This is equivalent to the well known 
removal of inhomogeneous broadening by the spin echo. 
Any pulse spin resonance sequence containing instanta- 
neous 7r/2 or TT-pulses can be described by Eq. (5). An- 
other important example is the class of Carr-Purcell se- 
quences used for coherence control (7r/2 — [r — tt — r — 

echo] repeat). 

III. DANGLING-BOND SPIN RELAXATION: 
DIRECT VS. CROSS-RELAXATION 

The presence of an inversion center in crystalline Si 
leads to weak spin-orbit coupling and extremely long spin 
relaxation times. The Ti for localized donor electrons in 
crystalline silicon can reach thousands of seconds at low 
temperatures.^^ This is in contrast to spin-lattice relax- 
ation of dangling-bonds in various forms of amorphous 
silicon where instead Tf*^ was found to range between 



one and a hundred milliseconds at the lowest tempera- 
tures (T = 0.3 — 4 K).^'^ The proposed theoretical expla- 
nation was that DB spin relaxation happens due to its 
coupling to phonon-induccd transitions of tunneling two 
level systems (TTLS) in the amorphous material. The 
TTLSs are thought to be structural rearrangements be- 
tween groups of atoms, that can be modeled by a double 
well potential [see Fig. 3(a)]. The TTLS assumption is 
able to explain several special properties of amorphous 
materials at low temperatures.^^ The DB spin couples 
to the TTLSs either through spin-orbit or hyperfine in- 
teraction, both of which are modulated by the TTLS 
transition. Note that the presence of a TTLS breaks the 
crystal inversion symmetry. 

We start by developing the theory of phonon-induced 
transitions for the TTLS,'^^ and the associated cross- 
relaxation of the DB spin. The Hamiltonian for a TTLS 
reads 

The energy scale e is a double well asymmetry, while 
A = AqC"^ is the tunneling matrix element between the 
states [A is related to the barrier height and its thick- 
ness, see Fig. 3(a)]. After diagonalizing Eq. (8) we ob- 
tain T^TTLS = diag{£;/2,-£;/2}, with E = V?TA2 
(for notational clarity we prime the Hamiltonians in the 
non-diagonal basis) . The coupling to phonons can be ob- 
tained by expanding the parameter e to first order in the 
phonon strain operator. 



^-E^^l^lKe-- + 4e-n, (9) 

leading to e ^ e + e'ii. Below we average over TTLS 
parameters with e S> A, so to be consistent we must 
assume the deformation parameter A' = 0. Applying this 
expansion to Eq. (8) and transforming to the diagonal 
basis we get 



"^TTLS-ph — 




Using Fermi's golden rule for dissipation into a phonon 
bath Tiph — '^gfy^qa\aq, we find that the transitions 
from +E/2 to —E/2 and vice- versa are given by 

r+ = aSA2 [nph(S) + 1] , (11) 
r_ - aEA^nph{E), (12) 

with phonon occupation number 

In Eqs. (11), (12) the parameter a depends on the ma- 
terial density p, sound velocity s, and deformation po- 
tential e' [a = (87r|e'p?i4ps5)-i]^ r^j^g ^^^^ Zeeman 



energy is denoted by Hbb = h-y^BS^^. To simplify 
the notation we define 5 = tvy^.B as tlie DB spin Zee- 
man energy. The coupling of the DB spin to the TTLS 
may be derived directly from the spin-orbit interaction 
Hso = aS^^ -{Exp), where S^^ is the DB spin op- 
erator, p is the DB orbital momentum, and E a local 
electric field. After averaging over the coordinate states, 
the resulting effective Hamiltonian becomes directly pro- 
portional to the magnetic field, a consequence of time 
reversal symmetry. For simplicity, we assume that E 
is perpendicular to the interface, and that the spin- 
orbit energy fluctuates by a certain amount A x 5 when 
the TTLS switches. This leads to the following effective 
Hamiltonian in the non-diagonal basis 

W'ttls-db = f (5?^ + ( ) > (14) 

where are raising and lowering operators for the DB 
spin. The dimensionless constant A will play the role of 
a small parameter in our theory. Transforming to the 
diagonal basis we get 

Wttls-db = § {Si- + S^-) ( t ) • ^''^ 

As a resuh of Eq. (15), the DB-TTLS cigenstates are 
admixtures between spin up and down. We may still label 
the eigenstates by their spin quantum number, provided 
we think of "f (J,) as having a large projection onto the 
pure spin up (down) state. The four level structure is 
shown in Fig. 3(b) in the hmit E ^ S and in Fig. 3(c) 
for E<^S. 

The total Hamiltonian is given by 

H = TixTLS+'^DB+'^TTLS-DB+'Wph+'WTTLS-ph- (16) 

Note that the first three contributions denote the dis- 
crete TTLS-DB states (a four- level system), the fourth 
is the energy bath (a continuum of phonon states) and 
the fifth is the coupling between the TTLS-DB to the 
phonon bath. The eigenstates of the first three contribu- 
tions may be calculated using perturbation theory, and 
the transition rates are straightforward to compute. The 
"direct" relaxation rate, corresponding to a DB spin-flip 
with the TTLS state unchanged is given by 

with [nph{S) + 1] ^ nph{S) for the reverse rate -D±x^±t- 
Note that Eq. (17) is proportional to A'*', reflecting 
the fact that a direct spin-flip may only occur together 
with a virtual transition to an excited orbital state. ^^'^'^ 
In our case this virtual transition is a "double-switch" 
of the TTLS, hence £> oc A* [terms independent of 
A in Eq. (17) cancel exactly. This general feature of 
a direct spin-flip process is referred to as "van Vleck 
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FIG. 3: (a) Effective double well potential for the tunneling 

two level system (TTLS). (b,c) Energy level structure for a 
dangling-bond spin (DB) coupled to a TTLS, for (b) _B > 5 
and (c) E <^ 6. (d) Coordinate system for the interaction of a 
dangling-bond located at rt with the donor spin, tpi denotes 
the angle formed by the donor-DB vector (dashed) and the 
external B field. 

cancellation" giving a simple explanation of why direct 
spin-flip rates are generally weak]. Moreover, Eq. (17) 

vanishes at i? = in accordance with time reversal sym- 
metry (the direct process couples a Kramers pair). 

The "cross" -relaxation rates, whereby the DB spin 
flips simultaneously with a TTLS switch are given by 

r_x = a\M+f{E + 6)np,,{E + 5), (18) 

= a\M+f{E + S)[np^(E + S) + l], (19) 

r_| = a|M_|'(S-5)nph(£;-<5), (20) 

F+i = a\M.\^{E~S)ME-S) + l], (21) 

where the sub-indexes label the level that the system is 
exiting, for example F+j = F+|^_|. Note that the final 
state is obtained from the initial state by changing the 
sign of the TTLS and fiipping the DB spin. The matrix 
element M± is defined by 

AfA 

M± = ^[\E±6\+6]. (22) 

Remarkably, this cross-relaxation process is not a tran- 
sition between Kramers conjugate states. As a result, the 
rates are qualitatively different from the direct process, 
particularly due to their magnetic field (S) and TTLS en- 
ergy {E) dependence. At low temperatures (fcsT -C 6), 
the direct rate always scales as £> oc 5^}"^'"^^ In con- 
trast, the cross-relaxation rate has two distinct behav- 
iors, depending whether E 6, or E <^ d. For E ^ 6, 
M± w AeA/E, and the F's are independent of magnetic 
field. For < (5 we get instead M± « 25AeA/E'^, and 
r (X (J'^ in contrast to the scaling of the direct rate. 
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Of extreme importance to our theory is to note that 
whenever the energy scales E and S are well separated, 
the direct rates are much smaller than the cross relax- 
ation rates. For £; » 5 we have D/T ~ {A/e)'^{S/E)^, 
while for E<^S we have D/F ~ {E/e)^{A/6)'^. The typi- 
cal assumption for amorphous semiconductors is A ^ e, 5 
and E « e.^° In this regime the direct rates are substan- 
tially weaker than the cross relaxation rates, except at 
the resonance point E = S. It is useful to list simple 
expressions for the cross-relaxation rates in the two most 
physically relevant regimes considered in this work. For 
low magnetic field S < ksT, E 6 hut with E/ksT 
arbitrary we have simply 



(23) 



Hence when the spin-orbit coupling parameter satisfies 
A <C f , the cross-relaxation spin- flips are much less fre- 
quent than the spin-preserving TTLS switching events. 
The opposite high magnetic field regime with E <^ 6 and 
E <^ ksT with d/ksT arbitrary leads to 

13 42 a2 

r+T«r_T^rT«4a— ^[nph(5) + l], (24) 

with the reverse rate Ti given by [nph((5) -I- 1] ^ fiph(^). 
Note that these T^^i rates are still much larger than the 
direct rates, since -D/r-fj^ ~ (A/(5)^. 

Finally, we discuss how the cross-relaxation rates are 
affected by the presence of phonon broadening in a non- 
crystalline material. In this case we generalize our theory 
by including a complex part to the phonon spectra, ojq = 
sq + i'^ph- The modified Eq. (22) becomes 



S2 + .,2 



^J{Elsf 



+ ■ 



{E±5f+^IJ2 



■7ph 



E±5\ 



iE±df+j, 



ph. 



(25) 



For amorphous Si we estimate 7ph ~ O.Qlsq.^^ For E 
S, 7ph ~ O.OIE may be comparable to S, and we see that 
M± is reduced by a factor of two, and an additional B 
field dependence results. 



IV. DANGLING-BOND SPIN RELAXATION: 
ENSEMBLE AVERAGE 

In order to evaluate the ensemble averages over TTLS 
parameters we must first determine the time-dependent 
correlation function for the four-level relaxation network 
described in Fig. 3(b,c). Using the notation of Eq. (3), 
the magnetic dipolar field produced by a single DB spin 
maps into a c-number fiz = 2hdipS^^ ^dipSi, with 
Sj = +1 (DB spin up) or = —1 (DB spin down). 
In the four- level system notation (+ J,,— J,.+ "f,— T) 
the vector x of dipolar fields assumes the values x = 



^dip(— 1, +1, + !)• In Appendix B we prove the con- 
venient identity 

S^t) = {[Vz{t) - fj,] [vM - fj,]) = X ■ p{t) ■ x^. (26) 

Here x^ = (xiwi, X2W2, ■ ■ ■), with Wi the equilibrium 
probabilities for the i-th level of the DB-I-TTLS network. 
The matrix p{t) — e~^* describes the occupation prob- 
ability for each level, and decays according to the relax- 
ation tensor A. Below we discuss the important analytic 
solutions for Sz{t) in the limit of small spin-orbit cou- 
pling, A < 1. 



A. Case £ > (5, (5 < ksT, E/ksT arbitrary 

In this regime the TTLS and cross-relaxation rates are 
simply related by Eq. (23). The time correlation function 
for the DB spin may be calculated exactly from Eq. (26), 
but for simplicity we show the result to lowest order in 

powers of A: 



Sz{t) 



hip 



*e-('^++''--^)*-h(l-*)e- 



rt 



(27) 



with a visibility loss given by 



^ = 



tanh^ {E/2kBT) 
cosh^ {E/2kBTy 



(28) 



and a thermalized DB spin relaxation rate given by 
2r- ^ . 2r+ ^ 



r = 



r+ -I- r. 
2aA^A^ 



-r+ + 



r+ + r- 
E 



sinh [E/kBTY 



(29) 



where we used E. Interestingly, Eq. (27) shows that 
DB spin relaxation happens in two stages: In the first 
stage the DB spin decays abruptly to a small visibility 
loss 5', with a rate set by the TTLS switch. During 
this first stage the TTLS levels ±E/2 achieve thermal 
equilibrium. In the second stage the DB spin relaxes 
fully with a much slower "thermalized" cross-relaxation 
rate T. For A <C 1 we may drop the oc A* contribution 
to Eq. (27). 

The theory developed above can be generalized to a 
single DB coupled to an ensemble of TTLSs, provided 
the TTLSs are not coupled to each other. In this case 
the rate Eqs. (29) and (37) are generalized to a sum of 
rates Fj relating to the i-th TTLS. Each exponential in 
Eq. (27) becomes ~ e^^ZiTit 'pj^jg happens whenever 
the DB-fTTLS network can be separated into discon- 
nected four level subspaces as in Fig. 3(b,c). 

We now proceed to average over disorder realizations 
of the amorphous material. We assume the following 
two-parameter distribution 



P{\,E) 



Pv 



E 



E„ 



(30) 
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for A e [0,A„,aJ, and E e [£^min, £^max]; P{\E) = 

otherwise. Note that the uniform distribution in A leads 
to a broad distribution of TTLS tunneling parameters 
A = Aoe""**. To our knowledge there are no estimates 
available for P, i?max, £^min close to an interface, only 
for bulk Si02. For the latter material the cnc;rgy den- 
sity of TTLSs per unit volume P has been estimated 
as P = 10^° — lO^^eV"^ cm~^, and typical values for 
the TTLS energy range arc ii'min/fcs ^ 0.1 K, and 
Emax/ks ~ 10 K.^^ Here we introduce a new parame- 
ter V with units of volume, denoting the effective range 
for TTLSs to couple to a each DB spin (for a Si02 layer 
of 10 nm wc estimate z; ~ 10'' nm"'). The exponent a is 
material dependent: While a « seems to be appropri- 
ate for bulk Si02,^^ it was found that bulk amorphous 
Si can be described by a = 0.1 — 0.4 or a = 1.2 — 1.5 
depending on sample preparation method (See Ref. 13 
and section IV C below). 

The average number of TTLSs coupled to each DB spin 
is given by 



J\f= J dxj dE F(A, E) 



PvEn 



1 



(31) 



This is also the number of thermally activated TTLSs 
at high temperatures, ksT » -Bmax- For lower tempera- 
tures satisfying £',„in <C fc^T <C -Emax, Eq. (31) is divided 
by cosh^ {E/2kBT), leading to 



ATt « PvE„ 



2kBT 



1+a 



(32) 



This is the number of thermally activated TTLSs inter- 
acting with each DB spin. For extremely low temper- 
atures ksT <C i?min this number will be exponentially 
small. 

Wc now turn to computations of the ensemble averaged 
DB spin relaxation rate, (F). At shorter times satisfying 
FMax* ^ 1, the DB spin magnetization {Sz{t)) decays 
linearly in timc.^^ The rate for this linear decay is equiv- 
alent to the (l/rf ^) rate measured for bulk amorphous 
silicon samples in Ref. 13. This is given by 

(f) = = /^^/ dEP{X,e)T{X,E). (33) 

At high temperatures ksT » E'max Eqs. (33) and (29) 
lead to 



(f ) = aA^AlksT 



An 



(34) 



The average DB spin relaxation scales linearly with tem- 
perature times the number of TTLSs surrounding the 
DB. 

At lower temperatures satisfying -Bmin ^ ksT <C -Bmax 
we have instead 



(r> = ?^(*BT)--/ 

J 



A R° 

/\max-^inax 
12 a27 



dx- 



sinh^ X 



= -T ^TT^^T, 

■^max-^ 



(35) 



showing that the DB spin relaxation rate will scale pro- 
portional to 

At the very lowest temperatures fc^r <C -Bmin there are 
no thermally activated TTLSs, therefore the mechanism 
of DB cross-relaxation is exponentially suppressed. Here 
other sources of DB spin relaxation may dominate [e.g. 
direct relaxation as in Eq. (17)], or the DB spin may 
not relax at all within the characteristic time scale of the 
experiment. 

Askew et al. measured average DB relaxation rates 
in bulk amorphous silicon at low temperatures (T = 

0.3 — 5 K).^'' Two different preparation methods, sili- 
con implanted with ^*Si, and silicon sputtered in a sub- 
strate, led to the experimental fit (F) cx T^.as ^-^q other 
preparation methods, silicon implanted with ^"^Nc, and 
silicon evaporated on a substrate led to T ' "' and T'^-^ fits 
respectively. Two different values of the magnetic field 
were studied (0.3 and 0.5 T), and no magnetic field de- 
pendence could be detected. The T and B dependence 
predicted by our model agrees with experiment provided 
a = 0.35 for the ^^Si implanted and sputtered samples, 
and a = 1.3, 1.5 for the ^°Ne implanted and the evapo- 
rated samples. It's perhaps expected that a is different 
for each of these because the density of TTLSs should 
depend on the way they were created. At high temper- 
atures, the linear in T behavior has been observed in 
amorphous silicon grown by evaporation.^'' 



B. Case E « 6, E « ksT, S/ksT arbitrary 



From Eq. (24) and Eq. (26) we get 



cosh^ {6/2kBT) 



For E <^ S the DB relaxation rate becomes 

F,+F,«4aA^|;^3,,,,^_^^. 

Its ensemble average is given by 



(36) 



(37) 



(rr+Fi) 



( E 



X I i±^]AA53coth 

1 — a 



IkeT 



, (38) 



where we assumed a < 1. For a > 1, the pref- 
actor in Eq. (38) is modified, but the scaling oc 
MS^ coth {5/2kBT) remains. 



C. Comparison to Ref. 13 

We now compare our results to the theoretical model 

proposed by Askew et al.}^ In their Eq. (5) the authors 
wrote the expression for F in the E ^ 6 regime using 
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free parameters D, Ad , C, N. In our work these are ex- 
plicitly related to microscopic parameters: D = e'e/E, 
M ^ -e'A/E, C = e/{2E), N = A/{2E). In Ref. 13 it 
is claimed that when the inequality ND/E ^ —CM/6 
is satisfied, the average DB relaxation rate scales as 



DB\-1 



rp2+ago (^^YiQ gp called Lyo and Orbach regime 



after Ref. 12). When this inequality is reversed, they 
obtained (rfB)-i oc T^+^J-^ (Kurtz and Stapleton 
regime, Ref. 11). Nevertheless, our result shows that 
these parameters are related by ND = —CM > 0, 
so this inequality is equivalent to S » E. Because 
Eqs. (7) and (8) of Ref. 13 are based on two conflicting 
approximations, 6 ^ E for the matrix element squared 
and S <^ E for the phonon density of states, their result 
needs to be corrected. We showed above that the av- 
erage DB relaxation scales instead as 5'^ coth.{S/2kBT) 
for (5 > and T'^+'^S° for S <s: E (the latter holds 
for i^min -C ksT <C £^max- For high temperatures 
ksT » £^max we get TJ"). The corrected results are 
in excellent qualitative agreement with the experimental 
data in Ref. 13. 

Ref. 13 assumes e = A = E/\/2 and averages E ac- 
cording to a density ^ E". This is in contrast to our 
averaging prescription that assumes instead A = Aoe""^, 
with Aq < emin and as a consequence e k, E. We as- 
sume A is uniformly distributed and the e density varies 
as ~ . This assumption is motivated by the wide distri- 
butions of TTLS relaxation rates observed in glasses, and 
is usually employed to explain charge and current noise 
in semiconductors. As we show below, the broader dis- 
tribution of DB relaxation times leads to 1// magnetic 
noise and non-exponential relaxation for an ensemble of 
DBs. 



V. MAGNETIC NOISE 

The total noise power for each DB spin is independent 
of the specific relaxation process and may be calculated 
exactly using elementary Boltzman statistics. The noise 
must satisfy the following sum rule: 



J — < 



(hi,) 



cosh-" {5/2kBT) 



(39) 



This shows that the noise spectrum is exponentially small 
in the high magnetic field regime 5 ^ ksT. For the oppo- 
site regime; 6 <C fc^T the total noise power is independent 
of temperature. However, as we show below, the spectral 
density ^^(w) may be temperature dependent when its 
upper frequency cut-off is temperature dependent. 



A. Case S » 5, <5 < ksT, E/ksT arbitrary 

In order to determine the noise spectrum, we must first 

extract the distribution of relaxation rates P{T) from 
Eqs. (29) and (30). Under the assumption that each DB 



spin is coupled to only one TTLS on average [i.e., ^ 1 
or At ~ 1, see Eqs. (31), (32)] we have 



p(r') = 



/ / dEP{\ E)6 (r(A. E) - F') 
JdXj dEP{\, E) ■ 



(40) 



Note that this is normalized to one according to 
/dr'P(r') = 1. It is straightforward to extend Eq. (40) 
to a larger number of TTLSs Ei, E2, ■ ■ ■, but the explicit 
calculation of P(F') becomes difficult. Below we will de- 
rive explicit results for the case of a DB spin coupled to 
a single TTLS on average. 

Using Eqs. (29), (30), and (40) we may evaluate the 
integral over A explicitly: 



P(F') = jj= j dEP{0,E) j dX 

dEP{0, E) 



6[X-Xo{E)] 

ldAU=^o(B) 



1 1 



2aA^E 



sinh {E/ksT) 



r' 



(41) 



Here Ao(£;) is the solution of r(Ao,-B) = V . The step 
function results from the fact that the delta function 
will "click" only when Aq {E) G [0, Amax] , or simply 
r' < 2a E / sinh {E/kBT). 



1. High temperature, ksT ^ -Emax 

In this case the theta function in Eq. (41) is always 

one for F' G [FminjFmax], with fmax = 2aA^kBT and 



min ^ 



^Fjnax- Therefore we have simply 



P(F') 



9 \ r' ' 



(42) 



for F' e [Fmin,F„iax], and P(F') = otherwise. As a 
check, note that / rfF'P(F') = 1 implies the relationship 

Amax = 5 In , as expected. 

The magnetic noise is given by 



S{^) = (hlJ / rfrp(r) 



F/tt 

,2 + r2 



4 Amax 



(43) 



for Fmin < uj < Finax. and S{uj) = for w > Fmax- For 
w < Fmin it saturates at 5(fmin)- Hence at the highest 
temperatures we have temperature independent magnetic 
1/f noise. 

The 1/f frequency dependence shows that the average 
magnetization of an ensemble of DB spins out of equilib- 
rium will decay non-exponentially with time t. At inter- 
mediate times satisfying f "5,^^ <C t "C f we may show 
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that the time correlation function (or cquivalcntly the 
ensemble average of the DB z-magnetization) satisfies^° 



1 - 



CE+ln(r^ax^) 



(44) 



(5?B(0)) 

This expression is valid after neglecting terms 
0(l/fmaxi)- Here Ce = 0.5772 is the Euler-Mascheroni 
constant. 



2. Intermediate temperature, Emin <C fesT <^ Err 



In this case Eq. (41) becomes 
1 \ + a f ksT 



p(r') 



2r Xr, 



E„ 



l+a 



The 



upper 



limit 



dx X 



of the 



'sinha; 



• (45) 



sinhx 



r'/r 



integral is determined from 
We solved this equation numerically, 



and showed that the result is well approximated by the 



analytic expression x^ 
approximation we get 



111(21^) 

\ ^-L max / 



Using this 



p(r') 



1 



2A, 



E 



l+a 



r' 



l+a 



(46) 

The distribution of relaxation rates has the same temper- 
ature dependence as the number of thermally activated 
TTLSs [see Eq. (32)], and possesses an interesting loga- 
rithmic correction with respect to the usual 1 /F' behav- 
ior. 

The logarithm correction in Eq. (46) increases the 
weight for smaller rates F', at the expense of decreas- 
ing the weight for higher rates. As a result the noise 
spectrum is better described by a 1//'' relation, with 
p > 1. Fig. 4 shows numerical calculations of S{ix>) 
for a = 0,0.35,1.5 (we assumed Fmin = 1 s~^, and 
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For a = 0, the noise is described 

1.5 



by a 1// fit, while for a = 1.5 a fit of 1// is more 
appropriate. Therefore at intermediate temperatures we 
have 



{hi. 



knT 



4'^max V-E'max 



l+a 



1 

\J\P 



(47) 



Note that X'^^.^ is determined from the normalization 
condition J (L)S{u)) = (/i^jp) for given Fmax/fmin- 



3. Extremely low temperature, ksT <C -Emin 

In this case f (A, E) is exponentially suppressed, and 
there will be no magnetic noise due to the DB+TTLS 
mechanism. If spin relaxation is dominated by the direct 
process [Eq. (17)], the noise spectra may still have the 
1 // dependence. Otherwise paramagnetic DBs may not 
contribute to magnetic noise at all. 




10"' 10° 10' 10' 10' lO'* 10^ 

CO (rad/s) 

FIG. 4: Magnetic noise at intermediate temperatures fcsT ^ 
S and Emir, < ksT < S^ax, for a = 0,0.35,1.5 (TTLS 
energy density exponents). The distribution of relaxation 
rates [Eq. (46)] contains a logarithmic correction, leading to 
S{cu) oc l/f, with p = 1.2 - 1.5. 



4. Calculation of (hip) 

Finally, we calculate the total noise power by averaging 
the DB distribution over the interface plane. We choose 
a coordinate system with origin at the interface immedi- 
ately above the donor spin. Define d as the donor depth, 
and Ti, (t>i the coordinates of the ith DB with respect to 
the interface [sec Fig. 3(d)]. The dipolar frequency shift 
produced by a DB spin aligned along the same direction 
as the donor spin is given by 



{hdip)i 



4(d2+r2)'/' 



(1 — 3cos^ Vj) 



(48) 



/idip is sensitive to the orientation of the external mag- 
netic field B = (sin6',0,cos^)S with respect to the in- 
terface. This enters through 



2 , id cos 6 + ri cos sin 6) 

cos = — 2 

d^ + rf 



(49) 



For = 0, the average over an uniform DB area 
density ctdb is given by 



(^dip) = <^DB / # / TdrH, 
Jo Jo 



10 Jo 
= 64^°^^- 



dipV) ' 



(50) 



VI. HAHN ECHO DECAY DUE TO 1// NOISE: 
COMPARISON TO EXPERIMENT 



The discussion above concluded that the following 
model for the noise spectrum is valid at high temper- 
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atures (fc^T > S and ksT > E^^^): 

{C'/Wmin < |a;| < Wmin 
C/\UJ\ Wmin < Iwl < W 
max 
Wmax < |w| < 00 

The prefactor C is given by 

(/idip) 37r CTDB 7e?l^ 



(51) 



c 



4A„ 



256 An 



rf4 



(52) 



Wc may calculate the Hahn echo response due to 1// 
noise using Eq. (5) with the filter function Eq. (7). If 
the inter-pulse time r is neither too long (so that Cmin = 

7rTCL)min/2 < 1) nor too short (so that Cmax = T^TCOmax/^ > 

1) we get 



,2 

min 



((T+(2r)) = expj - Cr^ 

X ( 3 - 4 cos (2Cinax) + COS (4Ci„ax) 



1 



4c2 



?3) 



after neglecting terms of order c^;^ and 1/c^ax- When 
Cmin ^ 0.1 and Cmax ^10 the echo envelope saturates 
and is well approximated by the simpler expression 



(ct+(2t)) « exp [-41n(2)CT2 



(54) 



that is independent of the low and high frequency 
plateaus assumed in Eq. (51).^*'^^ 

In the experiment of Ref. 3, each implanted Sb donor 
is a probe of magnetic noise from the interface. Because 
the implanted profile is inhomogeneous, the parameter C 
is different for each layer of donors a distance d below the 
interface. The experimental data was taken at S/ksT = 
0.3/5 = 0.06 « 1. Prom Eq. (52) we obtain 



(a+(2r))«e- 



,,, 6.25 nm 2 



(55) 



(56) 
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0.8 
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o 
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X ins) 

FIG. 5: Theoretical calculations (solid lines) and experimen- 
tal data (circles)"' for Hahn echo decay of Sb donors in the 
120 KeV implanted sample with the Si/Si02 surface. The 
theory is in reasonable agreement with the data when the 
theoretical parameter ^ ~ 0.2 [Eq. (57)]. The inset shows 
the Sb donor distribution measured by Secondary Ion Mass 
Spectroscopy (SIMS). 



400 KeV implant (annealed) 
SIMS 
Smoothed 




1000 2000 
X(|IS) 



4000 



FIG. 6: Same as Fig. 5 for the 400 KeV implanted sample, 
with the Si/Si02 surface. As pointed out in Ref. 3, the exper- 
imental data suffered from external field noise for r > 500 //s. 



C = (57) 

Amax 

In this approximation wc may fit the experimental data 
using a single dimensionless parameter ^, provided the 
distribution of Sb donors is well known. 

We used Eq. (55) together with the donor distribution 
measured by Secondary Ion Mass Spectroscopy (SIMS) to 
obtain theoretical estimates of Hahn echo decay relevant 
to the experiment of Ref. 3. Figs. 5 and 6 compares the 
theory with the 120 KeV and 400 KeV implanted sam- 
ples respectively, both with a Si/Si02 surface. A value 
of ^ 0.2 for the theoretical curves seems to be consis- 
tent with the experimental data. However, in the short 
time range the theoretical curve seems to decay slower 



than the experimental data, while at longer time inter- 
vals the theory seems to decay faster. This lack of agree- 
ment may be due to deviations from the measured SIMS 
distribution. The ultra-low donor densities were at the 
sensitivity threshold for the SIMS technique, hence the 
donor distribution is quite noisy [see insets of Figs. 5, 6 
- we used a numerically smoothed version of the SIMS- 
annealed data of Figs. 1(a), 1(b) of Ref. 3]. A higher 
probability density near the interface could in principle 
explain the faster decay at shorter times, while a deeper 
tail in the distribution could be responsible for the slower 
decay at longer times. 

The value for A^ax may be estimated from Amax = 
iln^Ea^") - llnfr^^ ~ 10. Combining this with 
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^ - 0.2 we get ctdb ^ lO^'^cm"^. 

We use Eq. (54) and the value ^ « 0.2 extracted from 
experiment to estimate the coherence time of a single 
donor located a distance d below the interface. This re- 
sults in 



T2{d) w 4 X IQ-^s 



(58) 



with T2 (d) inversely proportional to the square root of the 
DB area density. The 1// noise affecting a local magnetic 
probe a distance d from the interface is estimated as 



6.5 X 10"s-2 



10 nm 



and is directly proportional to (Tdb • 



VII. MAGNETIC FLUX NOISE IN SQUID 
DEVICES 



The SQUID (superconducting quantum interference 
device) is probably the most sensitive probe for mag- 
netism at the nanoscale. It consists of a superconducting 
loop interrupted by two insulating barriers (Josephson 
junctions). In this way it works as a magnetic flux-to- 
voltage transducer. SQUIDs are usually grown on top 
of a Si/Si02 substrate, therefore magnetic noise due to 
dangling-bonds within the SQUID loop will affect their 
performance as sensitive magnetometers. 

Our results on magnetic dipolar noise are easily trans- 
lated to flux noise in SQUIDs by substituting (/i^jp) for 
(^Totai) section V above. In order to get an order of 
magnitude estimate for (<I>xotai)' consider the flux pro- 
duced by a single magnetic dipole moment mo located 
at the center of a disk of radius R (the area enclosed 
by the SQUID loop). In Gaussian units this is given by 
$i = 27rmo/i?. Each dangling-bond contributes a dipole 
moment equal to mo = fi.7e/2 = he/ {2mi,c). Assuming 
an area density ctdb for the DBs leads to the following 
estimate for the mean flux squared: 



(*lotal) 



= 2.49 X 10"^^$^ [ctdb X (nm)2] , (60) 

where $0 = hc/2e is the flux quantum. The SQUID 
operates at very low magnetic fields (_B < 1 G), so the 
spin quantization direction is set by local inhomogeneities 
and is different for each DB. The angular average of spin 
quantization direction reduces Eq. (60) by a factor of 
3. Moreover, taking account of spins close to the su- 
perconducting wire and oriented along the SQUID plane 
increases Eq. (60) by ^ 3.^ As a result, Eq. (60) has the 
same order of magnitude as the calculation of Koch et al 
for loop sizes 10 — 500 ixm7 

At high temperatures (ksT 3> S and fc^T ^ i?max), 
the flux noise due to the presence of DBs in the plane 



enclosed by the SQUID loop is obtained by substituting 
(^dip) ^ (^Totai) in Eq. (43). The result is 



S,s,{uj) 



o-pB X (nm)^ 

-^max 



-12^0 



6.2x10-1^-^. (61) 



The value in brackets equals the parameter ^ used to flt 
our ESR experiment (Fig. 5). Using ^ w 0.2 we get an 
estimate for the flux noise contribution from an untreated 
Si/Si02 substrate: 



S^{u) « 1.2 X 10 



<I>2 
-12 ^0 



(62) 



Interestingly, this result has the same order of magnitude 
as the measurements of Ref . 6 using a small flux qubit as a 
probe of magnetic noise. A compilation of measurements 
of flux noise in SQUIDs was given recently by Ref. 7, 
where we see that S'(lIIz) lies in the range 0.1 — 100 x 
10~^^$o ^ wide variety of samples. 

Note that the high temperature condition ksT ^ S 
implies T ^ 0.1 mK for the low magnetic fields (^ 1 G) 
in SQUIDs. Unfortunately, there are no estimates of 
Emax for a Si/Si02 interface. For bulk Si02 the values 
-E'max ^ 10 K and -Emin ^ 0.1 K were estimated. We em- 
phasize that Eq. (62) is the maximum value for the noise, 
which saturates at ksT i^max- For ksT < Smax, 
Eq. (62) will be reduced by a factor {KBT/EynaxY^" , 
and the frequency dependence will change to 1/|w|p with 
p= 1.2- 1.5, see Eq. (47). 



VIII. NUCLEAR SPIN NOISE FROM A 
HYDROGEN PASSIVATED SURFACE 

Surface passivation with hydrofluoric acid drastically 
reduces the amount of dangling-bonds. Nevertheless this 
occurs at the expense of adding a large amount of hy- 
drogen nuclear spins. Here we investigate the magnetic 
noise mechanism arising due to the dipolar fluctuation 
of hydrogen nuclear spins at a perfect passivated Si-H 
surface. 

It is well established that spin decoherence of donors 
in bulk natural silicon is dominated by nuclear spin 
noise from the 4.67% ^^Si nuclear spins. The sam- 
ples studied here^ are known to have less than 0.1% 
of ^^Si isotopes, leading to a contribution of the order 
of ^ ~ 10 Hz. For a Si/Si02 sample, the fraction of 
oxygen isotopes with non-zero nuclear spin is even lower 
(0.038%), hence oxygen nuclear spins should be a minor 
contributor to magnetic noise at oxidized samples. 

We carry out a model calculation for the Si(100)H sur- 
face under the assumption that the hydrogen atoms are 
arranged in a canted-row dihydride phase with no orien- 
tation disorder, see Fig. 7.^" The truncated Hamiltonian 
for a single donor electron spin interacting with the hy- 
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FIG. 7: A Hydrogen terminated silicon surface is obtained 
after immersing an oxidized sample in a hydrofluoric acid so- 
lution. Here we show a Si(100)H surface with the hydrogen 
atoms forming a canted-row dihydride structure.^" The SiH2 
groups form a square lattice of side 5.43/\/2 = 3.84 A. 

drogen nuclear spin lattice at the surface is given by 

i i 

+ J2bij{Ii+Ij--'^IizIjz), (63) 

where li is the nuclear spin operator for the hydrogen 
atom located at position Ri with respect to the electron 
[jn = 2.66 X 10'* (sG)~* is the gyromagnetic ratio for the 
hydrogen nuclear spin]. Note that in Eq. (63) we have 
neglected the non-secular contribution of the electron- 
nuclear dipolar interaction. This approximation is valid 
only at higher external magnetic fields, jgB ^ vSi^i • 
For d = 10 nm, S > 0.1 T is necessary to satisfy this 
criteria. One can show numerically that the non-secular 
interactions produce a loss of visibility for the Hahn echo 
envelope scaling as ^ -D?/(7e-B)^.^* The electron- 
nuclear dipolar coupling is given by 

A = ^(l-3cos2^i), (64) 

where 6i is the angle between Rj and the direction of the 
external magnetic field. Each pair of hydrogen nuclear 
spins labeled hy i,j are mutually coupled by the dipolar 
interaction 

bij = -\^{^-^oosHij), (65) 

where 9ij is the angle between the B field and the vector 
Rij linking the two nuclear spins. 

The Hamiltonian Eq. (63) is directly mapped into the 
effective model Eq. (1) through the prescription rjz = 
Diliz- The noise spectrum [Eq. (3)] is then calculated 
using a "flip-flop" approximation, i.e. assuming a model 
that considers only flip-flop transitions between pairs of 
nuclear spins. In this approximation, the noise spectrum 
becomes^'' 

^M = ElS^['^('^ + ^^i)+'^('^-^^i)]' (66) 



with Ay = (D, - Dj)/4:, and E,, = 2^6?^. -FA?.. We 

dropped the inhomogeneous broadening term propor- 
tional to 5{oj) in Eq. (66) because it does not contribute 
to Hahn echo decay. 

Using Eq. (5) and Eq. (66) the Hahn echo envelope 
becomes 

(a+(2r)) = []e-4^-^^-^-'^'-'(-V^5+^), (67) 

i<j 

where sine a; = sinx/x. This result is identical to the 
lowest order cluster expansion derived in Ref. 20 through 
direct calculation of the spin echo response. Another 
way to derive Eq. (67) is to assume that the nuclear spin 
pair transitions are quasiparticle excitations with infinite 
lifetime. "^^ Eq. (67) is able to predict the Hahn echo de- 
cay without any phenomenological fitting parameter, in 
contrast to the traditional "Brownian motion" models 
developed previously. ■^^ 

Note that the magnetic noise due to nuclear spins is a 
linear combination of sharp peaks (delta fimctions), re- 
flecting the mesoscopic nature of the nuclear spin bath. 
Each delta function is a transition between discrete nu- 
clear spin energy levels. This is in contrast to the con- 
tinuous (Lorentzian) noise due to a single dangling-bond 
spin interacting with the phonon continuum. 

In order to plot a continuous noise spectrum we repre- 
sent the delta functions in Eq. (66) by normalized Gaus- 
sians with linewidth cr = 10^ s^*. Note that the Hahn 
echo decay is independent of the particular choice of a or 
the Gaussian lineshape provided r remains much smaller 
than 1 /a [in this case the Hahn echo envelope calculated 
by Eqs. (5), (7) with a coarse grained noise spectrum is 
very well approximated by the zero broadening expres- 
sion Eq. (67). 

Fig. 9 shows the nuclear spin noise spectrum from the 
point of view of a probe (a donor spin) lying 10 nm below 
the surface. Interestingly, we find that the noise spectrum 
is sensitive to the relative orientation of the external mag- 
netic field with respect to the surface. The noise has a 
global minimum for 9 50°. As shown in Fig. 10, this 
effect translates into a variation of about 50% in the elec- 
tron spin dccohercncc time T2 [T2 is obtained as the 1/e 
decay of the Hahn echo given by Eq. (67)]. This orien- 
tation dependence is surprisingly different than the one 
in bulk Si:P, see e.g. Fig. 8 of Ref. 14. Fig. 10 shows 
that T2 is minimized when = and maximized when 
9 « 50°, in contrast to bulk Si:P where precisely the op- 
posite was found. This special orientation dependence is 
the fingerprint of nuclear spin noise in a Si(100)H sur- 
face, allowing a clear identification of this mechanism in 
pulse spin resonance. 

Fig. 11 shows T2 as a function of the donor distance 
from the surface. Note that we find T2 > 10 ms for 
d ~ 30 nm, suggesting that this mechanism should not be 
playing a dominant role in the shallow implanted sample 
of Ref. 3 (120 KeV sample). For d > 100 nm, T2 is 
hundreds of milliseconds, so hydrogen nuclear spins are 
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FIG. 8: Magnetic noise spectrum for the oxidized interface as 
probed by a single donor spin a distance d below the interface. 
For comparison, we show the nuclear spin noise spectrum in 
a hydrogen passivated surface. 
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FIG. 9: Magnetic noise spectrum due to the Si(100)H surface 
as probed by a single donor spin 10 nm below the surface. We 
show the noise spectrum for four different angles 6, labeling 
the relative orientation of the external magnetic field with 
respect to the surface. 
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FIG. 10: Magnetic field angular dependence of T2 for a donor 

spin located 10 nm below a hydrogen terminated surface. 6 is 
the angle between the external magnetic field and the (100) 
direction. The resulting orientation dependence is quite dis- 
tinct from the one due to ^®Si nuclear spins in bulk natural 
silicon. "'34 
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FIG. 11: Spin decohcrcnco time T2 as a function of the donor 
distance from the surface, for a passivated Si(100)H surface 
containing no dangling-bonds (squares) and for a Si/Si02 in- 
terface containing a dangling-bond density equal to lO^'' cm~^ 
[Eq. (58)]. 



not affecting the 400 KcV implanted sample cither. The 
nuclear spin noise in a passivated surface may be further 
reduced by a factor of ~ 4 by using deuterium instead 
of hydrogen (the deuterium gyromagnetic ratio is 3.28 
times smaller than hydrogen). This results in donor T2's 
greater by a factor of two. 

For a perfect hydrogen passivated surface the theo- 
retical T2S are much longer than the values measured 
in Ref. 3. It is well known that c;hemical passivation 
of a Si(lOO) surface can not remove all dangling-bonds, 
in contrast to Si(lll) that usually removes nearly all 
dangling-bonds.^^ Therefore the dangling-bond mecha- 
nism might still be playing a role in the passivated sam- 
ples. Repeating the experiment for the Si(lll) surface 
could possibly yield even longer coherence times. The 
finite Sb density in these samples implies that the mu- 
tual interaction between donor spins (donor-donor dipo- 
lar coupling) might play a role, a mechanism of decoher- 
ence referred to as "instantaneous diffusion" .^^ We have 
confirmed this expectation by showing that the contribu- 
tion to T2 due to instantaneous diffusion is of the order 
of 0.3 and 1 millisecond for the 120 KeV and 400 KeV 
samples respectively. Therefore instantaneous diffusion 
might explain a fraction of the measured echo decay 
rates. Refs. 36,37 discusses a method for completely re- 
moving the instantaneous diffusion mechanism in a doped 
sample. 

With respect to SQUID devices, we remark that the 
nuclear spin flip-flop mechanism considered in this sec- 
tion does not contribute to magnetic flux noise (a flip- 
flop preserves the value of the magnetic moment for two 
nuclear spins, leaving the total flux unchanged). The sta- 
tistical fluctuation of individual hydrogen nuclear spins 
(due to a finite T^) should be extremely small because 
is usually hundreds of seconds or more. The nu- 
clear spin noise due to ensemble fluctuations of nuclear 
spins may be detected by SQUIDs under optimal circum- 
stances, see Ref. 38. 
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IX. DISCUSSION 

In summary, we developed a theory of magnetic noise 
due to spin-flips of paramagnetic centers at an amor- 
phous semiconductor-oxide interface. The mechanism of 
danghng-bond spin relaxation due to its interaction with 
tunneling-two-levcl systems and phonons of the amor- 
phous interface was discussed in detail. We also showed 
how these effects may be greatly reduced by surface pas- 
sivation with hydrogen. Substituting the paramagnetic 
dangling-bonds with a monolayer of hydrogen nuclear 
spins reduces the magnetic noise level by many orders 
of magnitude, as seen in Fig. 8. We related these results 
to decoherence of spin qubits in silicon as a function of 
their distance from the interface and flux noise in SQUID 
qubits. 

Our work generalizes and extends the model of 
dangling-bond spin-lattice relaxation in amorphous ma- 
terials originally proposed in Refs. 11-13. Particularly, 
we clarified the difi'erent temperature and magnetic field 
dependence as a function of the ratio between TTLS en- 
ergy splitting E and DB spin Zeeman energy 6. 

The theory of paramagnetic DB spin relaxation is sig- 
nificant for two recent proposals of single spin measure- 
ment based on spin-dependent recombination of conduc- 
tion electrons with dangling-bonds close to the Fermi 
level. In these experiments the time scale Tf*^ sets 
the limit on single spin measurement fidelity. To our 
knowledge there is yet no experimental study of Tf ^ 
at the Si/Si02 surface. We propose the measurement 
of magnetic field and temperature dependence of DB 
spin relaxation at short times [Eq. (35)] and the non- 
exponential decay at longer times [Eq. (44)] in order to 
validate our theoretical results and give a full character- 
ization of the free parameters. 

Our calculations provide benchmark values for the ul- 
timate coherence times of group V donor spin qubits im- 
planted in an actual device structure made from nuclear- 
spin-free silicon. Although the longest coherence times 
are in principle achievable with a perfect oxidized surface 
without dangling-bonds, the inevitable presence of a large 
density of these defects in real devices make surface pas- 
sivation an attractive alternative. Since each donor must 
be positioned close to an insulating interface in order to 
allow gate control of exchange,^ hyperfine couplings, 
as well as electron shuttling,'*^ the interface effects de- 
scribed here will play an important role in the material 
optimization of silicon devices exploiting spin coherence. 

Recently, ^^Si nuclear magnetic resonance experiments 
in polycrystalline silicon at room temperature were inter- 
preted using a model of magnetic 1/f noise. ^ The pro- 
posed mechanism was related to the charge fluctuation 
of trapping-centers at the surface of the microcrystals. 
Our work suggests that it is the spin-flip of paramag- 
netic DBs, not trapping-centers, that probably account 
for most of the 1/f noise observed in Ref. 4. 

Koch et al. proposed a model of 1// flux noise in 
SQUIDs based on electron hopping to localized defect 



sites, and concluded that a quite high trapping-center 
area density (5 x 10^^ cm~^) was required to explain flux 
noise in SQUID qubits.'' Our work suggests that the spin- 
flip of paramagnetic centers from the substrate may pro- 
vide an alternative explanation, based on a more physical 
paramagnetic dangling-bond density similar to the one es- 
timated in their work. 

We remark that a. C — V analysis of an unannealed 
Si/Si02 interface leads to an energy density equal to 
p' ~ 1013 eV-icm-2 (See Fig. 4 of Ref. 9). This im- 
plies that the SQUID substrate is contributing at most 
ksTp' ~ 10"'^'' cm~^ of trapping-center area density at 
T = 0.1 — 4 K. Nevertheless, the area density for param- 
agnetic DBs should correlate with U p' 10^^ cm~^. The 
value obtained here (10^^ cm^^) is a factor of 10 higher. 

STM experiments provide another way to estimate the 
trap energy density. In Ref. 44, a clean Si(lOO) surface 
was exposed to low pressure oxygen in order to pro- 
duce approximately a single oxygen monolayer. When 
the tip to surface voltage was ^ 1 V, ten to one hun- 
dred trapping-centers could be detected in a 65 X 65 A2 
region. This leads to an energy density in the range 
10^'' - 10" eV^cm-^, higher than the C - F measure- 
ments. 

The frequency and temperature dependence of flux 
noise in SQUIDs was measured a while ago in Ref. 5, 
using a wide variety of samples. These included sil- 
icon oxide substrates deliberately and not deliberately 
oxidized, as well as sapphire substrates. Some samples 
showed no temperature dependence, and the frequency 
dependence appeared to fit 1//^ , with p = 0.58 — 0.80. 
This frequency dependence can not be explained by our 
model. Nevertheless, the absence of temperature depen- 
dence may be explained by our model, provided the ma- 
jority of DBs are connected to one or more thermally acti- 
vated TTLS (or equivalently, E'max < ksT). An interest- 
ing question for future research is whether the interaction 
between DB spins can account for this discrepancy. 

This work establishes an important connection be- 
tween flux noise in SQUID devices and ESR studies of 
implanted donor impurities or dangling-bonds. As a re- 
sult, ESR characterization may play an important role in 
the pre-screening of novel materials for SQUID fabrica- 
tion. 
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APPENDIX B: CALCULATION OF THE DB SPIN 
TIME-DEPENDENT CORRELATION FUNCTION 



APPENDIX A: MAGNETIC 1// NOISE DUE TO 
DANGLING-BOND CHARGE 
TRAPPING-CENTERS 

Here we consider a mechanism of magnetic noise simi- 
lar to the McWhortcr model of current noise in semicon- 
ductor devices, that does not involve phonons.^° 

Dangling-bonds with energy close to the Fermi level 
act as charge trapping-centers, capturing electrons from 
interface states. Magnetic noise occurs because each time 
an interface electron tunnels into the DB, it produces an 
effective dipolar field in the donor spin, that is given by 
Eq. (48) divided by two. 

The tunneling rate for a trap located a distance z from 
the interface is assumed to be r(z) = r^axe"^/^. Here 
z measures the depth of the charge trap into the Si02 
dielectric, z = refers to a trap at the Si/Si02 in- 
terface. We have A = •^fWJJ^mFE^j'X ~ 1 A, where 
iJo ~ 4 eV is half the gap difference between the two ma- 
terials. Fmax ^ 10^ s~^ depends on the cross section for 
electron capturing. Assuming a uniform distribution for 
z in the interval [0, ^;max] leads to a distribution of rates 
equal to 



p(r) 



n 



I dr I 

I dz I 



(Al) 



for Fmin < F < Fmax, and zero otherwise. The noise 
spectrum is readily calculated as 



5(0;) 



{hl^) I drp(F)-i^, 



1 

4 - 



512 2„ 



(p'/esT) 



rf4 w' 



(A2) 



Therefore we have C/w noise with parameter C = 

^(^dip) I (S^max). In order to compare to Hahn echo decay 
data we define a parameter ^ similar to Eq. (57). This is 
given by 



A 



-p'A;BT(nm^). 



(A3) 



The area density for DB's with energy close to the Fermi 
level is estimated as yo'fcsT, with p' ^ lO^'' cm^^eV"^.^ 
The maximum possible value of ^ (at T = 5 K) is esti- 
mated from p'ksT < 10^° crn^^. Assuming Zmax ~ lOA 
we get ^ < 10~^. We remark that this maximum possible 
value for ^ is four orders of magnitude smaller than the 
value required to explain the experimental data of Ref 3 
(see Figs. 5 and 6). 



We are concerned with finite frequency fluctuations of 
the DB spin magnetic moment along the B field (z di- 
rection), therefore we will describe the dissipative kinet- 
ics of the DB-I-TTLS network considering only diagonal 
density matrix elements in the Bloch-Wangsness-Redfield 
theory. 

Define the propagator matrix Pij{t) as the diagonal 
density matrix element pait) subject to the initial condi- 
tion pim{0) = SijSmj- This is just the probability that 
the DB-I-TTLS will be at the state i at time t given 
that it was at state j at time t = 0. Note that here 
the indexes i labels one of the four DB+TTLS levels 
(+ L- i,+ t - T) [See Fig. 3(b)]. Furthermore, de- 



fine the matrix Pij{t) = Pij{t) 



Wi 



with Wi = Pii{oo) 



the equilibrium probabilities for level i. The steady state 
solution is then given by Pij{t oo) = 0. With this def- 
inition the rate equations for the 4x4 matrix p becomes 
simply 



dt 



p{t) = -K-p{t), 



(Bl) 



with initial condition Pij(O) = 5ij — Wi, and a relaxation 



tensor A defined as follows: For i ^ k, An- = —Tk^i, 
that's minus the rate for entering level i from level k. 
For i = k, An = J2j{^i) ^i^j^ t;he sum of rates for 
exiting level i. The case of our four-level system without 
direct transitions between Zeeman sub-levels is described 
by 



f r+ + F+x 
-r+ 


-r+i 



- + r_ 
-r 





-i 





r+ + F+i 
-r+ 



-r-1 



— r_ 

- + r 



-t/ 



(B2) 

The vector for equilibrium probabilities w = 

(w+l, w;_l) is the eigenvector of A with 

eigenvalue zero, satisfying X^jUij = 1. 

In the four- level system notation (-|- J., — J.,-|- |,— f) 
the vector x of dipolar fields assumes the values x = 
/idip(— 1) —1) +1) +!)■ The correlation function becomes 

Sz{t) = {[Vz{t)-Vz][vM-f].]) (B3) 
= ^Xi[Pij{t) -Wi]xjWj (B4) 

= X ■ p{t) ■ Xw, (B5) 

with Xu, = {xiWi,X2W2, ■ ■ .)■ Eq. (B5) together with its 
explicit solution p{t) = e"-^* allows exact calculations of 
the correlation function Sz{t). 
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